NASA Technical Memorandum 4752 


Biomechanically Induced and 
Controller Coupled Oscillations 
Experienced on the F-16XL 
Aircraft During Rolling 
Maneuvers 


John W. Smith and Terry Montgomery 
Dryden Flight Research Center 
Edwards, California 



National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 



ABSTRACT 

During rapid rolling maneuvers, the F-16 XL aircraft 
exhibits a 2.5 Hz lightly damped roll oscillation, per- 
ceived and described as “roll ratcheting.” This phenom- 
enon is common with fly-by-wire control systems, 
particularly when primary control is derived through a 
pedestal-mounted sidearm controller. Analytical stud- 
ies have been conducted to model the nature of the 
integrated control characteristics. The analytical results 
complement the flight observations. A three-degree-of- 
freedom linearized set of aerodynamic matrices was 
assembled to simulate the aircraft plant. The lateral- 
directional control system was modeled as a linear sys- 
tem. A combination of two second-order transfer func- 
tions was derived to couple the lateral acceleration 
feedthrough effect of the operator’s arm and controller 
to the roll stick force input. From the combined sys- 
tems, open-loop frequency responses and a time history 
were derived, describing and predicting an analogous 
in-flight situation. This report describes the primary 
control, aircraft angular rate, and position time 
responses of the F-16 XL-2 aircraft during subsonic 
and high-dynamic-pressure rolling maneuvers. The 
analytical description of the pilot’s arm and controller 
can be applied to other aircraft or simulations to assess 
roll ratcheting susceptibility. 

INTRODUCTION 

During rolling maneuvers and tracking tasks, many 
modem fighter aircraft tend to experience roll instabili- 
ties or damped oscillations. These oscillations are 
noticeable at initiation and recovery following large 
primary roll control commands. Large roll rate com- 
mands require high gradient gains and high roll accel- 
eration to achieve and to recover from the roll in a 
desired manner. The periodic motions that often follow 
are sustained by the pilot/vehicle coupled interactions. 
Pilots and investigators generally refer to these interac- 
tions as either pilot-induced oscillations (PIOs) or “roll 
ratcheting.” The brief discussion that follows attempts 
to distinguish between the two. 

Usually, PIOs predominate during tracking tasks or 
following large, rapid control inputs. For example, dur- 
ing a tracking task, the pilot is actively minimizing the 
error or command signal. If the product of the pilot’s 
gain and the system gain increases sufficiently, then the 
pilot/vehicle closed-loop system tends toward neutral 
stability. Neutral stability PIOs most often develop at a 
crossover frequency within the range of 0.5 to 1 .0 Hz. 


Roll ratcheting predominates during roll recovery 
and is generally associated with or caused by an inter- 
action of the dynamic characteristics of the 
neuromuscular limb system, the aircraft, and the 
mechanical controller. The biomechanical nature of the 
phenomenon is similar to a “bobweight” in its 
feedthrough effect and produces oscillations within a 
frequency range of 2 to 3 Hz. The magnitude and 
damping are governed by the arm weight, the controller 
weight, the wrist interface, and the control system gain 
(primarily the commanded roll gradient gain). 

During the F-16 XL-2 functional test flight phase, a 
minimal number of rapid rolling maneuvers and recov- 
ery tests were conducted. At the end of the final func- 
tional test flight, the copilot, maneuvering and 
maintaining control from the aft seat, initiated and 
completed three rapid 360° rolling maneuvers at Mach 
0.9 and an altitude of 1 1 ,600 ft. During the roll recov- 
ery portion, as the copilot attempted to capture a 
desired roll angle, the aircraft exhibited a considerable 
amount of lightly damped roll oscillations. These oscil- 
lations continued through more than 10 cycles at a fre- 
quency of 2.5 Hz. Particularly important, the roll 
controller force transducer (measured at the controller) 
recorded the same characteristic oscillations as other 
aircraft control and response parameters. Based on the 
flight experience and confirmed feedthrough effects 
existing at the controller, the observed phenomenon 
was believed to offer a good example of roll ratcheting 
for analysis and reporting. 

Similar phenomena were encountered during other 
flight test programs. 1-5 The Air Force conducted F-16 
XL rolling maneuver evaluations. At Mach 0.9 at an 
altitude of 5200 ft, their pilots experienced roll ratchet- 
ing. Tracking studies using a variable-stability aircraft 
showed that roll ratcheting tendencies increase as aug- 
mented roll damping is increased and that pilot rating 
increases correspondingly. 4, 5 Using simulation and 
analytical studies where tracking tasks were generated 
by comparing position or target feedbacks to a forcing 
function, the problem has been defined and investi- 
gated extensively. 6-10 

In addition to presenting the actual flight data, this 
report also analyzes the biomechanical nature of the 
“bobweight” feedthrough effect and associated phe- 
nomenon. Linear models of the aircraft dynamics, con- 
trol system, and the arm, wrist, and controller were 
combined to model the phenomenon. 

A mathematical model and description for the pilot’s 
arm and controller were developed. Extensive use was 
made of existing data to define the wrist and arm 



dynamic characteristics. 10 These characteristics 
include wrist damping, spring constants, pilot’s arm 
equivalent weight, frequency, and damping. Previous 
tests had been conducted with the pilot partially 
restrained on a vibrating platform while performing 
single-axis tracking tasks. 10 Various controllers were 
firmly mounted to the vibrating platform. From these 
tests, the biodynamic constants and coefficients for the 
arm and wrist interface compliance were determined. 
The combination of limb dynamics, controller dynam- 
ics, and their equivalent weight effects are referred to 
as the manipulation dynamics. The mathematical 
description of the manipulation dynamics was concate- 
nated with and integrated into the aircraft aerodynam- 
ics and augmentation systems. 

A linear mathematical model for the integrated 
aircraft aerodynamics and augmentation systems was 
derived. The aircraft derivatives used to define the plant 
were available from previously reported flight 
results. 3,11-14 Wind-tunnel data for individually mea- 
sured control derivatives were also available. 14 The 
complete set of dimensional derivatives, as defined in 
appendixes A and B, was assembled as matrices in a 
state-space format. Using a FORTRAN control pro- 
gram, 15 open-loop frequency responses and a typical 
time history were calculated. 

This report describes the F-16 XL-2 roll ratcheting 
in-flight experience. Mathematical and theoretical 
analyses are presented to describe and predict the 


characteristics of the roll ratcheting phenomenon 
caused by a “bobweight” feedthrough effect. The flight 
test results and analytical predictions reported here 
were conducted at the NASA Dryden Flight Research 
Center, Edwards, California. 

VEHICLE DESCRIPTION 

The following sections provide a general description 
of the basic F-16 XL aircraft and aircraft control 
systems. Instrumentation and recording systems used 
to obtain the flight data were developed by NASA 
Dryden. 

Aircraft 

Two F-16 XL aircraft were designed and built by 
General Dynamics (Fort Worth, Texas) for a U.S. Air 
Force contract. The aircraft were initially considered 
prototypes, emphasizing the multirole fighter capabil- 
ity. At the start of the flight test phase, these aircraft 
were formally referred to as demonstrators. Both air- 
craft are modifications of the F-16 A/B aircraft. The 
forward and midfuselage have been extended by 56 in. 
The design planform incorporated a cranked arrow 
wing. The number 2 aircraft, of concern in this report, 
is a two-seat fighter using an F-110 GE-129 engine 
(General Electric, Lynn, Massachusetts) at the time of 
flight 12. Figure 1 shows the relative size and location 



Figure 1 . The F-16 XL-2 aircraft surfaces and pertinent control units. 
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of the primary and secondary control surfaces and rele- 
vant control system units. Table 1 shows the pertinent 
physical characteristics and geometry of the aircraft. 

Flight Control Systems 

The F-16 XL aircraft derives total vehicle control 
from aerodynamic surfaces. The surface positions 


depend on the primary control input and the feedback 
response as measured by the aircraft rate sensors. 

Aerodynamic Surfaces 

Flightpath control, in both pitch and roll, is achieved 
by a blend of the aileron and elevon surfaces. The sur- 
faces are faired along the trailing-edge segments of the 


Table 1. The F-16 XL-2 physical and geometric characteristics (ref. 3). 


Wing 
Area, ft 2 
Span (theo.), in. 

Aspect ratio 
Taper ratio 
Leading-edge sweep 
Inboard 
Outboard 
Incidence 
Dihedral 

Mean aerodynamic chord, ( c ) ft 
Elevon area, ft 2 
Aileron area, ft 2 
Leading edge flap Area, ft 2 
Vertical tail 
Area, ft 
Span, in. 

Aspect ratio 
Taper ratio 
Leading-edge sweep 
Bicon airfoil, percent 
Root 
Tip 

Rudder area, ft z 

Speed brake area, ft 2 
Control-surface authority (normal to HL) 
Elevons 
Ailerons 

Leading-edge flaps 

Rudder 

Speed brakes 

Control-surface rates, deg/sec at zero load 
Elevons 
Ailerons 

Leading-edge flaps 
Rudder 

Speed brakes, per panel 
Pertinent physical locations eg ~ 45% 
Lateral acceleration sensor 
Aft cockpit controller 


663.259 

388.84 

1.583 

0.128 

70° 

50° 

-0.65° at BL 41.50, -4.104° at 136.1 and tip 
0 ° 

24.7 

44.055 

29.445 

18.197 

54.75 

101.00 

I. 294 
0.437 
44.5° 

5.3 

3.0 

II. 65 
14.26 

30° up, and down 
20° up, 30° down 
6.4° up, 36.5° down 
30° left and right 
60° up and down 

60 

60 

31 

90 

30 

FS 324 
FS 118 
FS 162 
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wing. Directional control is accomplished by a conven- 
tional rudder mounted on the trailing edge of a single 
vertical stabilizer. 

Fly-By- Wire and Primary Control 

A full fly-by-wire analog system is used in all three 
axes. The system is fully quadruple-redundant. Electri- 
cal signals from the analog flight control computers go 
directly to the integrated servoactuators. These actua- 
tors provide the force and position for the primary con- 
trol surfaces. The electromechanical actuators are 
single-fail-operate components. 

Secondary Systems 

Leading-edge flaps (LEFs) tire located outboard on 
the cranked portion of the wing. The LEFs are sched- 
uled to operate symmetrically as a function of angle of 
attack and Mach number. Limited LEF differential trim 


input is available in conjunction with the normal aile- 
ron roll trim: 

ALEF (maximum) = 3.0° 

0 ( 1 ) 
for q c > 2176 lb/ft 2 

The LEFs also provide partial spin recovery by oper- 
ating differentially (ALEF = 40°) for angles of attack 
greater than 35°. 

Speed brakes are located on the aft portion of the 
fuselage inboard of the elevons and in the horizontal 
plane. The brakes operate by rotating into the vertical 
directions to provide deceleration and additional spin 
recovery capability. 

Sidearm Controller and Augmentation 

The F-16 XL airplane has the same two-axes 
minimum displacement sidearm controller used on pro- 
duction F-16 aircraft. Figure 2(a) shows the relative 



(a) Location of the pedestal-mounted controller within the cockpit. 
Figure 2. Pitch and roll minimum displacement sidearm controller. 
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cockpit location and the size and shape of the control- 
ler. Figure 2(b) shows the wrist and elbow rests pro- 
vided on F-16 aircraft. Use of the rests is optional; the 
rests can be stowed. Most pilots use only the elbow rest 
during rolling maneuvers. 

By applying force to the controller, the pilot initiates 
pitch and roll inputs to a rate command augmentation 


system rather than to controlling surface positions 
directly. Longitudinal force input commands a blend of 
pitch rate and normal acceleration. The pitch feedback 
loops attempt to provide response characteristics 
invariant with flight conditions. Lateral force inputs 
command roll rate according to a roll command gradi- 
ent schedule (fig. 3). The gradient consists of linear 












constant rate gains/lb over the ranges shown with a 
maximum commanded roll rate of 300 deg/sec. A force 
input of 9 lb and larger produces a maximum gain of 
30 (deg/sec)/lb. The force is actually determined by a 
transducer measuring displacement, Fy c = Fy c (x c ). 
While force and displacement are well-correlated stati- 
cally, the phase lag between applied force and resulting 
displacement is important in any transient analysis. A 
breakout force of ±1 lb is required for any transducer 
output, after which the transducer output is essentially 
linear up to a maximum displacement of 0.25 in. at 
20 lb (a spring constant of 80 lb/in or 960 lb/ft). Total 
surface position command is proportional to the error 
signal, e, as shown by the functional roll control system 
diagram (fig. 4). Favorable roll-yaw coordination is 
provided through a scheduled and function-generated 
aileron-rudder interconnect. 

A typical rate damper stability augmentation system 
is implemented in the yaw axis to accomplish both 
dynamic and static stability. Aircraft dynamic stability 
is augmented by a washout yaw rate feedback. Static 
directional stability is augmented by lateral accelera- 
tion feedback. 

A more complete and detailed description of the 
flight control, its design and final configuration philos- 
ophy, and the surface rates and authorities has previ- 
ously been given. 13, 16 


Instrumentation 

The F-16 XL-2 instrumentation was designed, 
installed, and maintained by the NASA Dryden Flight 
Research Center. Most of the data originate from flight 
control sensors and are obtained from the avionics data 
bus. The data presented in this report were recorded 
from pulse-code-modulation telemetry. The aircraft 
rates, surface positions, and aircraft system functions 
were recorded at 50 samples/ sec. Tables 1 and 2 show 
flight control accelerometer locations. 

FLIGHT TEST EXPERIENCE 

Many modem fighter aircraft experience roll 
oscillations. The occurrences follow rapid rolling 
maneuvers, particularly those performed at high- 
dynamic pressures. 

Roll Ratcheting 

During F-16 XL-2 flight 12, which in part was a 
functional check flight, the copilot performed three 
fairly rapid 360° rolling maneuvers from the aft cock- 
pit. The copilot reported severe roll ratcheting particu- 
larly at high commanded roll rates and in the recovery. 



Figure 4. Basic F-16 XL function roll control system for the flight conditions listed in table 2. 
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Table 2 shows the flight conditions typical of all three 
maneuvers. 


Table 2. Flight conditions, mass properties, and 
dimensions. 


Mach number 

0.9 

Altitude, ft 

11,600 

Dynamic pressure, lb/ft 2 

790 

True airspeed, ft/sec 

964 

Angle of attack, deg 

Angle of attack and pitch attitude, 

1.77 <a< 4.0° 

deg 

3 

Center of gravity, % c 

45 

Weight, lb 

28,318 

Impact pressure, lb/ft 2 

957.3 

Impact pressure/static pressure 
Vehicle moments of inertia, slug ft 2 

0.7 

i x 

22,352 

l y 

104,027 

I z 

123,293 

7 « 

Longitudinal distance, ft 

(A/C accelerometer location, 

-559 

relative to eg and body axis) 
Vertical distance, ft 

(A/C accelerometer location, 

17.0 

relative to eg and body axis) 
Controller location 

0.0 

Longitudinal distance, ft 

13.5 

Vertical distance, ft 

1.0 


Figures 5 and 6 show time histories of the aircraft 
response and control parameters. Figure 5 shows the 
time response of the lateral controller force, bank 
angle, rolling angular rate, and lateral acceleration at 
the sensor for the three rolling maneuvers and recovery. 
The direction of the roll is indicated. Initially, the copi- 
lot banked right to approximately 60° and rapidly com- 
pleted a full 360° roll to the left. The copilot continued 
the roll to the left at a slow rate and stabilized in a left- 
bank attitude for approximately 5 sec. From a wings- 
level attitude, the copilot completed another 360° left 
roll and stabilized the aircraft near 0° for approxi- 
mately 10 sec before initiating a full 360° roll to the 
right. 


At the recovery portion of all three rolling maneu- 
vers, where |F Kc (jc c )| > 9 lb, the aircraft exhibited a 
considerable amount of roll ratcheting. The data show 
the oscillations are generally superimposed on a 
commanded and derived roll response. Consequently, 
for one-half cycle of each oscillation, the motion 
appears to the pilot as a periodic hesitation in the roll 
response. Thus, the terminology “roll ratcheting” is 
very descriptive of the phenomenon displayed. 

Pilot Remarks and Comments 

The copilot who performed the rolling maneuvers 
has a considerable amount of time and experience in 
flying F-16 aircraft. From a discussion with the copilot 
following the flight and a review of the flight data, the 
subsequent impressions were obtained, and the 
copilot’s remarks and comments are summarized as 
follows. 

The copilot was definitely aware of the ratcheting 
response. The copilot’s primary attention was focused 
on the nose of the aircraft, so impressions were formed 
from the cyclic behavior and roll rate hesitation. The 
cockpit motion was of little concern during the occur- 
rence. The copilot also thought that the body, torso, and 
shoulders were restrained sufficiently to be minimally 
responsive. The elbow rest was in use while the rolls 
were executed. The copilot felt the forces were low or 
unchanging at the arm during the ratcheting accelera- 
tion experienced. This observation and comment is 
understandable because the controller force transducer 
was responding to inertia forces of the arm and 
controller. The calculation of the controller displace- 
ment resulting from the observed variation in controller 
force supports this explanation. The minimum dis- 
placement controller has a stiffness defined by 

F v (x r ) = 20 lb; x = 0.25 in. (2) 

l c c 

The measured controller force variations (from fig- 
ure 5(a)) are 

A Fy = 2.5 lb peak-to-peak (3) 

Then Ax c = (2.5/20) * 0.25 = 0.031 in. peak-to-peak as 
indicated in the Aircraft and Controller Cycle 
Response section. This small controller displacement 
seems to explain the lack of pilot observation. 
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(c) Roll rate. 



First 20 sec — 






Aircraft and Controller Cyclic Response 

By comparing the data from the three rolls (fig. 5), it 
is evident that the most extensive ratcheting occurred 
during recovery of the first rolling maneuver 
(10<f <20sec). The discussion and analyses that 
follow consider the first rolling maneuver and recovery 
as distinctive representations of the displayed 
phenomenon. 


Figure 6 shows time histories of the first 20 sec of the 
data shown in figure 5. These data present a complete 
360° left rolling maneuver and recovery. All the aircraft 
response data relevant to the lateral control and the 
appropriate control surface positions are shown. The 
most interesting time span of the time histories exhibit- 
ing the roll ratcheting is from 13 to 17 sec. These data 
show the existence of 10 cycles of oscillatory motion 
lightly damped and at a frequency of 16 rad/sec or 



(a) Stick force (from the stick position transducer). 



(b) Bank angle. 

Figure 6. The F-16 XL-2 expanded time histories during rolling maneuvers at flight conditions listed in table 2. 
10 



approximately 2.5 Hz. This frequency falls within the 
range of frequencies considered typical of roll ratchet- 
ing as reported in various references. 1-3 

Note that the lateral controller force, which is a func- 
tion of controller transducer output calibrated in lb, is 
very active during the roll ratcheting recovery 
(fig. 6(a)). This trace portrays a very lightly damped 
system oscillation and suggests that the controller is 
actually being excited through a “bobweight” effect 


that is sustained by lateral acceleration. At the initiation 
of the recovery, the force level |Fy c | was greater than 
9 lb. Figure 3 shows that the controller gradient dou- 
bles at a force level of 9 lb. The phenomenon observed 
is easily excited and sustained at high force levels. 
Fully saturated roll inputs ( | Fy c >17 lb) would not 
provide a gradient gain for the feedthrough process. 
However, during a typical recovery from a saturated 


P. 

deg/sec 




(d) Lateral acceleration (from the accelerometer in the front cockpit). 
Figure 6. Continued. 


11 



condition, the roll ratcheting oscillations are likely to 
occur. 

The time history data indicate a reduction in the 
oscillatory amplitudes as the controller force level 
decreases (fig. 6). The estimated peak-to-peak varia- 
tions of the distinctive aircraft responses during the roll 
ratcheting segment are as follows: 

Frequency range = 2.24 Hz < / < 2.83 Hz 
= 14 < co < 18 rad/sec 
co = 16 rad/sec (average) 

Aa y ~ 0.09 g 
Ap~2\ deg/sec 
A Fy (x) m 2.5 lb 

1 c L 

Ax c = 0.031 in. 

Ap - Apco = 336 deg/sec 2 

A dy(p) ~ 0.18g (at the controller) 

The peak-to-peak roll response, Ap = 21 deg/sec, 
resulting from a controller force variation of 
AFy c (.x c ) = 2.5 lb shows an amplitude of 8.4 (deg/sec)/ 


lb for a frequency of 2.5 Hz. This change is approxi- 
mately 28 percent of the commanded steady-state gain 
response (30 (deg/sec)/lb) for a controller force level of 
9 lb and greater. The acceleration variation measured at 
the sensor indicates a maximum of Aa Y = 0.09 g. This 
variation gradually decreases after the first four cycles. 
Little, if any, response is evident in the yaw rate, r, time 
histories. All primary control surfaces exhibit the ratch- 
eting oscillation. The rudder is particularly active and 
indicates a peak-to-peak displacement of more than 4° 
as a result of a y feedback. 

Although the accelerations are vector quantities, a 
rough estimate of the pilot’s arm weight can be 
obtained by adding the peak-to-peak magnitudes of the 
two components. This approach assumes that the arm 
can be modeled as a point mass at the controller 
(termed an equivalent weight). 

Aciy = Aa y + Ap — ■£=-?', 1 7 “ 1ft 
Y c r s g 57.3 z 

= 0.09 + 0.18 
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(e) Yaw rate. 
Figure 6. Continued. 






(h) Leading-edge flap position. 



(i) Rudder position. 
Figure 6. Concluded. 


because 


A Fy = 2.5 lb (5) 

1 C 


then 


AFy 

W T (e) = - = 9.3 lb (6) 

1 A a y 


In modeling the pilot’s arm, wrist, and hand and the 
sidearm controller, rotational as well as translational 


elements exist. The model converts the rotational ele- 
ments into translational elements by defining 
equivalent weights that include moment of inertia 
terms. The equivalent weight of the arm, W @ (e), is 
equal to the total equivalent weight of the arm and con- 
troller, W T (e ) , minus the weight of the controller, W c . 
The controller weighs 1 .25 lb. 

W @ (e) = W T (e)-W c = 9.30-1.25 

(7) 

W @ (e) - 8.05 lb 
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ANALYSIS AND DISCUSSION 

Previous roll ratcheting investigations and analyses 
considered the pilot not only as a connective link but as 
an active element of the closed-loop process. The fol- 
lowing discussion considers the pilot’s role as the 
active controlling element. However, the major investi- 
gation and analysis focuses on the force feedthrough 
effect resulting from lateral acceleration and the par- 
tially restrained controller and arm masses. 

Roll Ratcheting Investigation, Scope, and 
Definitions 

Studies of roll ratcheting oscillations encountered by 
high performance aircraft have considered the pilot’s 
role in a closed-loop process. 4, ^ The biomechanical 
feedthrough effect caused by the product of the lateral 
acceleration and the sum of the arm and controller 
masses have also been studied. 10 

Figure 7 shows two simplified illustrations depicting 
the control process in deriving a roll rate, capturing a 
roll attitude, or maintaining a bank angle. Figure 7(a) 
shows a simulated forcing function. The pilot, while 
continuously observing the desired error signal, 


actively tries to minimize the signal. How well the task 
is accomplished becomes a qualitative assessment and 
may be quantized as indicated for a comparison with 
another test configuration. The variables in most cases 
are the pilot or operator, the controller, and the con- 
trolled element. Analytically, the operator is modeled 
by a transfer function composed of a gain constant, a 
time delay, and some form of equalization. The con- 
troller is represented by its dynamic response charac- 
teristics. The controlled element is defined by the 
aircraft transfer functions. 

Figure 7(a) depicts a task system of an in-flight 
investigation using a T-33 variable stability aircraft 4, 5 
exemplifying the study of the pilot’s role. During track- 
ing tasks, the investigators varied such parameters as 
time delays, prefilter lag-lead models, and augmented 
roll mode damping. One reported finding was that roll 
ratcheting tendencies increase with augmented roll 
damping and the subsequent oscillations resulted in 
poor pilot ratings. This finding is logically consistent 
with the conclusions of this report. Increasing the aug- 
mented roll damping, decreases steady-state roll rate 
response to a pilot input. Large pilot inputs or a high 
gain on the pilot input are, therefore, required to mini- 
mize the error signal, which is a function of roll rate 



C(p, 4)) 


950178 


(a) System with a forcing function or task reference; operator controlling to minimize an error signal. 



Bobweight, feedthrough loop 


C(p, <t>) 


950179 


(b) System to achieve desired roll rate and roll attitude. 

Figure 7. Basic human operator system concepts for acquiring a prescribed task. 
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and bank angle. The large inputs result in large initial 
roll and lateral acceleration responses, which lead to an 
increase in roll ratcheting oscillations. Large inputs 
also result in increased stick gain (via nonlinear stick 
gearing). Such gains also lead to an increase in roll 
ratcheting susceptibility. 

The present investigation simply considers the pilot 
performing in-flight rolling maneuvers with no minimi- 
zation task required. Figure 7(b) shows this control 
process. First, the pilot applies a roll control force to 
the controller, Fy p , to achieve the desired roll rate. 
Then, the force is decreased or reversed near the 
desired roll attitude. Rolling the aircraft generates tan- 
gential acceleration components at the controller com- 
posed of, l x r, l z p and aircraft lateral acceleration, a y . 
The product of the equivalent weight and sum of the 
accelerations causes an additional force input, A F, that 
is superimposed on the reference force input, 
( Fy + A F). This reaction results in a closed-loop 
feedthrough effect that produces an overall reduction in 
gain margin. This closed-loop process is generally 
referred to as a “bobweight” effect. 

Assuming the arm is modeled as a cylinder with uni- 
form mass distribution rotated about an endpoint, then 
the equivalent weight of the arm at the controller is 
equal to one-half the weight of the arm, W @ , plus the 
moment of inertia effect. The moment of inertia effect 
is Wfe/3. 


W 

W 9 (e)m ^®+0.33W@ (8) 

Using the equivalent weight of the arm as estimated 
in equation 7 from the flight data, then 


have previously been given. 17 

The arm component 

weights for a 165-lb adult male cadaver, expressed as a 
percentage of total body weight, are as follows: 

Entire arm 

6.59 percent 

Upper arm 

3.44 percent 

Forearm and hand 

3.15 percent 

Forearm 

2.26 percent 

Hand 

0.89 percent 

Based on the F-16 XL pilot’s weight, 160 lb, the 
comparable arm components weights are estimated to 

be 


Entire arm 

10.54 lb 

Upper arm 

5.50 lb 

Forearm and hand 

5.04 lb 

Forearm 

3.62 lb 

Hand 

1.42 lb 

Assuming that the elbow rest was ineffective during 


the oscillations, then the equivalent weight at the con- 
troller would be determined from the weight of the 
entire arm ( W@= 10.54) as follows: 


W 

W @ (e) = -^ + 0.33 W @ 

10 54 

= — — + 0.33*(10.54) 

2 

W @ (e) = 8.75 lb 

Subsequent analytical calculations vary the equivalent 
weight of the entire arm arbitrarily between the two 
extremes (3.93 < W @ (e) $ 10.75) . 


8.05 = 0.83 W @ (9) 

W@= 9.7 lb 

If the elbow rest completely immobilized the upper 
arm, then the equivalent weight would be determined 
from the forearm and hand. The flight-determined arm 
weight of 9.7 lb appears to be approximately twice the 
expected weight for forearm and hand (approximately 
5 lb). Apparently, the elbow rest was only partially 
effective in restricting the motions of the upper arm. 

The total body weights of four cadavers and the 
weights of body components including arm segments 


Integrated Modeling Description 

Figure 4 shows a functional block diagram of the 
F-16 XL roll control system. Pilot input, Fy , is indi- 
cated as the reference input to the controller. A force 
input initiates the roll command through the force and 
roll gradient and is summed with the augmentation roll 
rate feedback. The result generates a roll response 
through the plant aerodynamics that is proportional to 
the error signal. Figure 4 shows the aileron-to-rudder 
interconnect and the angle of attack and roll rate (a p 
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term) configurations on the F-16 XL airplane. The 
input/output relationships for the plant aerodynamics 
are described in appendix B. Lateral-tangential accel- 
erations at the controller are generated because of the 
roll response and the additional lateral degrees of free- 
dom. These accelerations are computed in g. Along 
with the respective equivalent weights of the controller 
and pilot’s arm, these accelerations provide a force- 
feedthrough process back into and through the control- 
ler. The dynamic response of the controller was deter- 
mined from release tests with various weights added. 
The results are presented in appendix C. The natural 
frequency of the controller is well-separated from the 
ratcheting frequency. Consequently, the controller 
dynamics have a small effect on the ratcheting 
response. The response of the controller at the ratchet- 
ing frequency can be determined by the second-order 
response curves found in any dynamics text. 
For oo/co n < 0.3 and £ = 0.05, the gain is approxi- 
mately 1.0. 

The manipulation dynamics are shown in the biome- 
chanical feedthrough model (fig. 8). This particular 
illustration was adapted from a description previously 
given. 10 An extensive amount of experimental and 
analytical research regarding the biomechanical forces 


derived from accelerations was conducted. 10 This 
report relies heavily on those efforts and measure- 
ments. The complete formulation of the second-order 
systems and equation characteristics are derived in 
appendix D. 

Figure 9 shows a mathematical model of the 
biomechanical nature of an operator’s arm and wrist 
and the controller. The spring constant, damping coeffi- 
cient, and time constants assigned to the respective 
second-order system equations were derived from the 
vibration platform and experimental tests previously 
performed. 10 By letting s — » 0, the steady-state gains 
become 


X. 

i 






0.72 

0.28 


(ID 


or the total steady-state force to the roll command 
gradient is 

lim F y = 0.72 [F Y -W m (e)a Y )-W a y (12) 
s-»0 Y c ^ p Y ' J c 


Airframe 



Stick or 
controller 


Grip 


Limb 


Body 

semirestra'ined 
to the airframe 

950180 


Figure 8. Feedthrough mechanical controller model with lateral force and acceleration input. 
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Figure 9. Biomechanical model and mathematical description of the operator’s arm and controller dynamics 
(appendix D). The controller spring gradient, K c . . . , cancels out. 


If the operator’s hand is removed from the controller, 
then F Y , K t , and D ( = 0. Expectantly, this leaves the 
controller weight, w c , as the only force feedthrough 
loop into the aircraft roll command system. However 
because F v = 0, the roll command gradient gain is 

1 P 

very small (fig. 3). Consequently, little, if any, 
feedthrough effect would excite the forward loop. 

The equivalent weights (fig. 9) are considered 
representative and within the normal range of a typical 
arm and controller. The equivalent controller weight, 
, is considered equal to the actual weight of the con- 
troller at the disconnect location. The equivalent 
weights are used instead of mass and inertia because 
the calculated or measured tangential acceleration at 
the controller are conveniently expressed in g units. 

Mathematical Computations 

The flight control system (fig. 4) was linearized about 
the same flight test conditions (table 2) as the 
aerodynamics for the plant. A unique nonlinear closed- 
loop system, called a dual lag, is implemented on the 


F-16 XL airplane. This system provides a different 
command response depending on initial input, control 
bias, and frequency. Consequently, the control bias 
must be considered when determining the input/output 
relationship. Appendix E shows the development of an 
averaged linear transfer function for a sinusoidal input 
at the flight-experienced roll rate condition. 

CO _ 10(3 + 40) 13 

Rl (3 + 6)(3 + 67) 

The open- and closed-loop total system characteris- 
tics were determined by using a FORTRAN control 
program. 1 ^ Figure 10 shows an open-loop frequency 
response of the tangential acceleration/command force, 
a Y x /Fy . For the observed in-flight ratcheting fre- 
quency, 16 rad/sec, the computed system shows a gain 
margin of 27 dB and a phase margin of 70°. At the 
crossover frequency, co co = 26 rad/sec, the gain margin 
is 35 dB. Note that a gain increase of 10 dB exists near 
the ratcheting frequency. This gain increase undoubt- 
edly increases the tendency of the system toward insta- 
bility when coupled with the manipulation dynamics. 
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Figure 10. Open-loop tangential acceleration frequency response (conditions: § C /F y = 30 (deg/sec)/lb. No 
manipulation dynamics or time delay at the flight conditions shown in table 2). 
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Figure 1 1 shows a comparison of the vector 
components of the tangential accelerations at the con- 
troller using the same open-loop computation (except 
(0 = 15.6). The tangential accelerations are caused by 
roll acceleration, p, and the lateral-translational 
acceleration, a y + l x r. A ay resulting from p lags the 
translational acceleration by approximately 180°, and 
the amplitude of the p vector is 2.17 times greater than 
the translational vector ,f( a y< f). The resultant vector 
is also shown and has a phase angle of -107° with a 
magnitude of 0.047 g. 

Biodynamic Response 


Figure 12 shows the open-loop frequency responses 
for four equivalent total weight conditions: W T (e) = 
5.18, 7.25, 10.00, and 12.001b. For W T (e ) =5.18 lb, 
the open-loop system indicates a sufficient gain margin 
of 15 dB at (a co = 20 rad/sec. Progressively increasing 
the equivalent total weight decreases the crossover fre- 
quency and decreases the gain margin. At an equivalent 
total weight of 12 lb, the system becomes very lightly 
damped at a frequency of approximately 10 rad/ sec. A 
nominal and perhaps typical value for an equivalent 
total weight is believed to be approximately 10.0 lb, 
which means that the equivalent arm weight about the 
shoulder is as follows: 


Figure 8 shows the arm, wrist interface, and control- 
ler feedthrough model. Manipulation dynamic charac- 
teristics are derived and developed in appendix D. 
Figure 9 shows the numerical values. The total system 
is thus described by inserting the manipulation transfer 
functions into the aircraft and functional control system 
diagram (fig. 4). Using this completely integrated 
model and applying the mathematical control program 
connective procedures, 15 the dynamic characteristics 
and time responses of the total system were obtained. 


W @ (e) = W T (e)-W c = 8.75 lb (14) 

The equivalent weight includes a moment of inertia 
effect; 

W & (e) = (entire arm)/2 

(15) 

+ 0.33 W @ (moment of interia) = 0.8 lb 



Figure 11. Open-loop vector components of the resultant tangential acceleration at the controller (§ C /Fy = 30 
(deg/sec )/lb; to = 15.6 rad/sec. No manipulation dynamics or time delay at the flight conditions shown in table 2). 
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Figure 12. Open-loop frequency response as a function of total equivalent weight (fig. 9) (conditions: manipulation 
dynamics with no time delay at the flight conditions shown in table 2). 


or the actual weight of the entire arm is 

w@ = 10.5 1b (16) 

An arm weight of 10.5 lb is believed to approximate 
the pilot’s arm weight and yet be conservative enough 
to provide sufficient damping. 

Figure 13 shows a computed time history for a 1-lb 
step input at F y ■ Eigenvalues show the mode to have 


a natural frequency of 14.3 rad/sec and a damping ratio 
of 0.15. This frequency is within the range observed 
from the flight data. To determine the effects of time 
delays, 50 msec of delay was added in the forward loop 
of the system (fig. 4). This delay changed the natural 
frequency and damping ratio to co^ =13.2 rad/ sec and 

C = o.io. 


21 




A«Y t (P)- 

g 




Time, sec 


950185 


Figure 13. Analytically obtained time history in response to a unit step input (conditions: W T {e ) = 101b, 
<j ) c /F y = 30 (deg/sec)/lb, no time delay, at the flight conditions of table 2). 
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CONCLUSIONS 

The F-16 XL-2 airplane is a modified version of the 
production F-16 aircraft. Modifications include extend- 
ing the forward and midfuselage by 56 in. and chang- 
ing the wing to an advanced cranked arrow type. In 
place of horizontal tails, the F-16 XL-2 airplane derives 
total pitch and roll control by elevons and ailerons 
mounted along the trailing edge and throughout the 
span of the wing. The aircraft has a fly-by-wire analog 
flight control system approximately identical to the 
early F-16 aircraft. A minimum-displacement, two- 
axes, sidearm controller is used to command the pitch 
and roll rates. 

During a functional check flight, the copilot, flying in 
the aft seat, completed three rapid rolling maneuvers. 
Upon recovery from each maneuver, the copilot 
reported roll ratcheting oscillations that lasted through 
10 cycles at a frequency of 2.5 Hz. The roll command 
gradient gain was at its maximum of 30 (deg/sec)/lb 


initially and remained near the maximum throughout 
the oscillations. 

The following conclusions were reached from flight 
test results and the analytical studies regarding the 
phenomenon: 

1. During the roll oscillations, the controller trans- 
ducer showed oscillations of the same extent, 
phase, and continuation as other aircraft roll 
response parameters. The data support the hypoth- 
esis that the oscillations were commanded through 
the controller. 

2. Based on the controller force data, the maximum 
peak-to-peak controller displacement during the 
roll ratcheting oscillations was very small, Ax c < 
0.032 in. This displacement indicates small arm 
motions. The copilot observed no arm forces or 
motions and believed that there were none. 

3. All primary control surfaces followed the roll 
ratcheting oscillations induced by tangential accel- 
erations at the controller. 
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4. The open-loop, computed tangential response of 
the basic aircraft, a y/ F Y c > shows about 10 dB of 
gain increase near the roll ratcheting frequency. 
Coupling with the manipulation dynamics causes 
the closed-loop system to be prone to ratcheting 
instability. 

5. The tangential component of the acceleration at 
the controller caused by rolling acceleration, l z p, 
is approximately 2.2 times larger than the compo- 
nent caused by aircraft lateral accelera- 
tion, dy+ l x r. This observation indicates that roll 
ratcheting depends largely on p and on the dis- 
tance of the controller above the roll axis. 

6. The dynamic characteristics of the controller have 
only a slight effect on roll ratcheting. The ratchet- 
ing frequency is well- separated from the natural 
frequency of the controller determined in 
appendix C. 

7. Both crossover frequency and gain margin 
decrease when equivalent arm weight, W @ {e), 
increases. 

8. Based on estimates from one in-flight rolling 
maneuver, the wrist and elbow rests appear to be 


ineffective in restraining the arm from the small 
motions associated with roll ratcheting. The esti- 
mate of equivalent arm weight based on flight data 
more nearly matches the weight of the whole arm 
than just the forearm and hand. 

9. Following a step input of Fy p = 1 lb with a roll 
command gradient gain of 30 (deg/sec)/lb and an 
entire arm weight of 10.5 lb, the integrated system 
exhibits a slightly damped oscillation: £ = 0.15 at 
a natural frequency of 14.3 rad/ sec. 

10. Adding a time delay of 0.05 sec in the forward 
loop decreased the damping and frequency to C, = 
0.10 and co n = 13.2 rad/sec. 

1 1 . A technique for predicting roll ratcheting has been 
developed. The method requires an augmented air- 
craft model, the derived aircraft accelerations at 
the controller, and the pilot arm weight. The 
method and equations can be applied to other 
aircraft and included in simulation studies. 

Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, October 28, 1994 
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APPENDIX A 

NOMENCLATURE 

Abbreviations 

ARI 

aileron-to-rudder interconnect 

FS 

fuselage station, in. 

GE 

General Electric 

ISA 

integrated servo actuator 

LEF 

leading-edge flap 

PIO 

pilot-induced oscillation 

TF 

transfer function 

XL 

extra-long 

Symbols 

A 

matrix coefficient, ( A ' = C _1 A) 

a Y 

lateral acceleration at the center of 
gravity, g 

a Y c 

lateral acceleration at the controller, g 

C 

°y 

lateral acceleration at the sensor, g 

s 

a Y T 

tangential acceleration at the controller, j 

B 

matrix coefficient, ( B ' = C -1 B) 

b 

reference span, ft 

C 

matrix coefficient 


rolling-moment coefficient, rolling 
moment/qsb 

C l 

roll-damping derivative, 

p 

C/d(pb/(2V)), rad" 1 

C l 

rolling moment due to yaw rate, 

l r 

dC/d(rb/2V), rad” 1 

c i 

'p 

effective dihedral derivative. 

dC^/dp, deg -1 

c i 

\ 

differential aileron roll effectiveness 

derivative, 3^/38^, deg -1 

c i 

l & 

e 

differential elevon roll effectiveness 

derivative, 9Cy/98 f , deg -1 

c i 

‘6 

r 

rolling moment due to rudder deflection, 

dC/d8 r , deg -1 

c n 

yawing-moment coefficient, yawing 
moment/qsb 





c 

eg 

D 

D c 

D t 



/ 

I 


X ’ 


xz 


Y' l Z 


yawing moment due to roll rate, 
dC n /d(pb/2V), rad- 1 
yawing damping derivative, 
dC n /d(rb/2V), rad" 1 
directional stability derivative, 
dC n /d P, deg -1 

yawing moment due to aileron deflection, 
ac„/a8 a> deg- 1 

yawing moment due to elevon deflection, 
dC n /d8 e , deg -1 
rudder effectiveness derivative, 

ac„/a5 r , deg -1 

generalized control system output 
side force coefficient, side force/qs 

side force due to roll rate, 
dC Y /dp, deg- 1 

side force due to yaw rate, 
dCy/dr, deg -1 

side force due to sideslip, 

ac r /ap, deg -1 

side force due to aileron deflection, 
dC y /d8 a , deg -1 

side force due to elevon deflection, 
dC Y /dS e , deg" 1 

side force due to rudder deflection, 
dCy/d 5 r , deg -1 
mean aerodynamic chord, ft 
center of gravity, % c and FS 
damping, lb sec/ft 
damping of the arm, lb sec/ft 
damping of the controller, lb sec/ft 
damping of the wrist and hand interface, 
lb sec/ft 
equivalent 
force, lb 

lateral controller force, lb 
lateral controller force, lb 
lateral controller force pilot input, lb 
roll-to-rudder gain schedule, J{M) 
function of parameter inside parentheses 
vehicle moments of inertia, slug ft 2 
product of inertia, slug ft 2 
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" @ 

W @ (e) 

W„ 


imaginary part of Laplace transform 
variable, rad/sec 

gain constant 
arm spring constant, lb/ft 
controller string constant, lb/ft 
wrist spring constant, lb/ft 
roll rate gain 
longitudinal distance, ft 
vertical distance, ft 
aircraft mass, (lb-sec 2 )/ft 
mass of arm, (lb-sec 2 )/ft 
mass of controller, (lb-sec 2 )/ft 
rolling angular rate, deg/sec 
roll acceleration, deg/sec/sec 
static atmospheric pressure, lb/ft 2 
dynamic pressure, lb/ft 2 
impact pressure, lb/ft 2 
generalized reference input 
yawing angular rate, deg/sec 
reference planform area, ft 2 
Laplace transform variable, (c + jc o), 
rad/sec 
time, sec 

true airspeed, ft/sec 
weight, lb 

weight of the arm, lb 
equivalent weight of the arm, lb 
weight of the controller, lb 
total weight of the controller plus added 


A @ arm transfer function denominator 

A controller transfer function denominator 

c 

8 differential aileron deflection, 

<VV /2 ' deg 

8 differential elevon deflection, 

( 5 e R - 8 ^ )/2 ’ deg 
8 r rudder deflection, deg 

0 pitch altitude, deg 

a real part of Laplace transform, rad/sec 

<|) bank angle, deg 

§ C /Fy roll rate command gradient, (deg/sec )/lb 

A$c/A.Fy roll rate gradient gain, (deg/sec)/lb 

to frequency, rad/sec 

(O co crossover frequency, rad/sec 

oo ^ damped frequency, rad/sec 

00 n natural frequency, rad/sec 

co r ratcheting frequency, rad/sec 

II absolute value 

II dB amplitude ratio, dB 

Sign Conventions 

Trailing edge down is positive for 8 a , 8 fl , 8 e , 
& e . Trailing edge left is positive for S r . Accelerations 
are positive forward, right, and down. Angles and 
rates are positive right wing down, nose up, and nose 
right. Angle of attack is positive with the nose (longi- 
tudinal axis) above the velocity vector. Angle of side- 
slip is positive with the nose left of the velocity vector. 


W T (e) 

weight, lb 

total equivalent weight of the arm and 

Subscripts 

X 

controller, lb 
output quantity 

@ 

arm 

X@ 

arm displacement, ft 

C 

controller 

x c 

transducer displacement, in. and ft 

D 

damped 

a 

angle of attack, deg 

i 

wrist interface 

P 

angle of sideslip, deg 

L 

left 

8 

control surface deflection, deg 

n 

natural 

e 

error signal 

o 

trim condition 

; 

damping ratio 

P 

pilot 

A 

small change 

R 

right 
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T total 

X, Y, Z vehicle body axis 

x tangential 

Dot over quantity denotes first derivative with respect 
to time. 

Dimensional Coefficients 

The dimensional coefficients are defined as follows: 
L - ^ C 

b ~ , x S 

n b - 

p I z 
Y = f 

r P 

qSb^ r 

P 2 VI Y l P 


N = 

P 2VI Z n P 

Y d = + sina ■ 

p 2m V 2 Y p 

L = ^ C 

S a .S f .8 r l x ‘s a .s e ,s r 

N *«*"*r = q r^ Cn K\K 

Y K>KK = £/ Cy kkk 


_ qSbl 
' - 2 V//', 

AT = 

r 2V/ Z "r 

y = qSb c Y -cosa 
r 2m V 2 *> 


sina 


-cosa 
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APPENDIX B 

DERIVATIVES AND AIRCRAFT PLANT 
FORMULATION 


C„ = -0.0006 

"*r 

Cv = 0.001 


In general, the aircraft lateral-directional aerody- 
namic derivatives were obtained from either an average 
or the mean value of available flight data for a Mach 
number near 0.9. 3, 9 ’ 10, 1 1 Wind-tunnel data were used 
in obtaining the control effectiveness of the elevon and 
aileron control derivatives. 12 

The coefficients listed in equation (B-l) were 
obtained with reference to the wind or stability axis. 
These coefficients provide the basic aerodynamics for 
generating the plant output and the integrated aug- 
mented aircraft response. 

C, = -0.0016 

'p 

C = 0.002 

"P 

Cy = -0.012 

y P 

C, = -0.24 

p 

C„ = -0.02 


C / r = 

c_ = 


c 


c, = 


(y p = sin a) 
0.10 
- 0.2 

(F r = -cosa) 
-0.0013 

-0.0005 

0.002 

-0.0011 

-0.0007 

0.002 

0.0003 


(B-l) 


First, it was necessary to translate the coefficients 
from the stability to the body axis and then convert to 
the dimensional coefficients (eq. (B-2)) using the flight 
data and the physical characteristics (table 2). 


_ 4 sb r 

- 77S 


No 


N r 


- 4 Sb r 

h "P 

- 3Lr 

- mV L r p 

2 VI x 
_ qSbl r 

2 VI z n P 

= sin a 


L 

N 

Y 

L 


_ qSb 

” 

q Sb 


KKK i x c h a , 5,.8 r 


8 a - S r- 5 r I z n S a , 5 f . 5 r 

- Mr 

5 a > 6 e> 5 r mV K 8 a ,S f ,8 r 

_ q sb2 r 

- wr x c >, 


) Y. 


qSb-~ 


C n r 

Z r 


= -cos a 


(B-2) 


Using the state-space notation, the lateral-directional 
equations can be arranged initially in a matrix format 
as 
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The complete computational procedure, including 
the axis rotation, is accomplished by using a FOR- 
TRAN program developed at the NASA Dryden Flight 
Research Center. Numerical values for the matrix coef- 
ficients A ' , B ' , H, and G at the example flight condi- 
tions are presented as equation (B-5). 


-.4172 .0523 -.9986 .0333 

-73.172 -3.0602 1.2365 0. 

15.1645 -.0369 -.4511 0. 

0. 1.0 .0524 0. 


0.071 .0348 .0695 

-54.285 14.267 -45.26 
-4.152 -4.588 -5.659 
0. 0. 0. 


(B-5) 


1. 0. 0. 0. 

0 . 1 . 0 . 0 . 

H = 0 . 0 . 1 . 0 . 

0. 0. 0. 1. 

0. -1.567 29.92 -.998 


0. 0. 0. 0. 

0 . 0 . 0 . 0 . 

G = 0. 0. 0. 0. 

0. 0. 0. 0. 

29.93 0. .53 0. 


The final linearized numerical matrix coefficients 
(eq. (B-5)) are a suitable aerodynamic state-space 
representation for integration into the FORTRAN con- 
trol program. 13 Output-dependent variables of the 
aerodynamic plant provide the rates and accelerations 
necessary for augmenting the aircraft. 


28 



APPENDIX C 

MINIMUM DISPLACEMENT 
CONTROLLER DYNAMICS 

The F-16 XL-2 aircraft was constructed with a dual 
cockpit during the modification and assembly. Both 
crew members are provided with identical two-axes, 
pedestal-mounted, minimum displacement controllers 
that are located on the right side of the cockpit. The 
controller is positioned slightly forward of the center of 
gravity of the pilot’s body to permit a semi-relaxed arm 
position. 

In both axes, the controller senses displacement mea- 
sured through the transducer displaced by force up to a 
maximum displacement of ±0.25 in. The transducer is 
calibrated in lb. A breakout force of 1 lb is required 
before any displacement results. The spring constant up 
to the maximum deflection is approximately 960 lb/ft 
(fig-C-1). 



Figure C-l. Spring gradient of the minimum displace- 
ment controller. 

To get some idea of the dynamic response of the con- 
troller as a function of mass, extra weights were added 
by wrapping lead strips around the controller grip. The 
controller, W , weighs 1.25 lb. 

The controller was pulled to maximum deflection 
and released. Figure C-2 shows a typical response for a 


total weight, of 5.18 lb (W added = 3.93 lb). For 
this particular weight and test, the controller exhibited 
a damped frequency, (O d , of 76.6 rad/sec at a damping 
ratio, £ , of 0.09 following the release. 



Figure C-2. Typical controller response following a 
Fy = 20 lb release: Wj =5.18 lb, period = 0.082 sec, 
co d = 76.6 rad/sec, £ = 0.09 and K c = 960 Ib/ft. 


Figure C-3 shows a summary of this test and two 
additional test weight conditions. For the controller 
alone (W = 1.25 lb), the release results in no over- 
shoot. A dashed line is drawn through the damping 
ratio data points. A constant value of damping force 
(D = 2.24 (lb-sec)/ft was used to estimate the damp- 
ing ratio from 5.18 to 9 lb. The solid line on the fre- 
quency figure indicates the natural frequency by using 
co n = jK c g/W T . The two data points were the 
damped frequency obtained from the response for the 
indicated added total weight. 
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Figure C-3. Controller frequency and damping as a function of total weight, W c + W(added). 




APPENDIX D 

BIOMECHANICAL MODELING AND 
MANIPULATION DYNAMICS 

The following description, derivations, and modeling 
aspects of the “bobweight” effects rely heavily on the 
experimental and analytical studies. 10 A large vibration 
platform was used on which a control operator was 
restrained at the seat to the platform. The platform was 
subjected to lateral accelerations ranging in frequencies 
from 0 to 10 Hz. The operator performed isolated-axis- 
displayed control tasks. Two types of center-stick con- 
trollers were mounted on the platform: a spring stick 
and a stiff stick. These controllers were investigated 
throughout the lateral control tests. From these experi- 
ments, the most significant modal responses were 
derived from the biomechanical nature of the arm/ 
muscle coupled through the hand/grip interface to the 
controller. 

Figure D-l shows a biomechanical model, showing 
the lateral forces and accelerations, and a mathematical 
description of the dynamic elements. 10 The directions 
of the forces, accelerations, and displacements are 
defined as positive to the right. Consider a minimum 
displacement controller, K c = 960 lb/ft, with aircraft 
accelerations measured in g and the mass properties 


defined in terms of equivalent or effective weights 
(fig. D-l). The forces and elements can be summed and 
resolved as two equations with two unknowns, x c and 

x @- 


[m c s 2 + (D ( + D c )s + K l + K c ]x c 

(D-l) 

~(D i s + K i )x @ = W c (-a Y ) 

- (D'S + K i )x c 

+ [m@s 2 + ( D i + D @ )s + K i + £@]x@ (D-2) 

= W @ (e)(-a Y ) + F Yp 


or 


[A c + (D'S + K i )]x c - (D-5 + K t )x @ 

= W c(-“Y ) 

- (D ( s + K i )x c + [A @ + (D ( .s + K-)]* @ 
= W @ (e)(-a Y ) + F Yp 


(D-3) 


(D-4) 


where 


Airframe 



Stick or 
controller 


Grip 


Limb 


Body 
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Figure D-l. Feedthrough mechanical controller model with lateral force and acceleration input. 10 
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A c = m c s 2 + D c s + K c 
A@ = m @ s 2 + D @ s + K @ 

m c = W /S 

m @ = W @ (e)/g 

Solving equations (D-3) and (D-4) for x @ and com- 
bining gives the following equation for the biodynam- 
ics: 


For simplification, the initial and end blocks can be 

AFy 

thought of as unity where K c = ■ - - (fig. (C- 1 )). 

The feedback block, D t s + K i /(D i s + K { ) + A @ , 
can be further moved to the forward loop and simpli- 
fied to a block diagram with F v and a Y as the inputs 
and Fy as the output or input to the roll command 
gradieiU(fig. (D-3)). 


(D-5) 


F Yp (D iS + + (D,j + KfKW^e) + W c ) + A @ W c 

(DjS + AT-)(A @ + A c ) + A @ A C 


However, A c is much larger than A @ ; therefore, 
A@ + A c ~ A c . For the lower frequency range of inter- 
est, co < 30 rad/sec, A c can be approximated by the 
controller spring constant, K c . Equation (D-6) then 
becomes (for a positive a Y ) as follows: 

A * T 



960083 


Figure D-3. Equivalent representation of figure (D-2) 
with pilot input. 


F Y p ( D i s + K i) /a Y x 
“ ^ c [(D t .5 + ^ ( ) + A @ ] 

(D., + ^.)(W @ ( g) + lV c ) 

K'UDiS + KJ + Aq) 

K c [(D i s + K i ) + A @ ] 

Equation (D-6) can be represented as a feedthrough 
system with tangential acceleration a y at the control- 
ler as the input and Fy as the output (fig. D-2). 


Measured values for the coefficients and gains were 
determined from vibration tests as (ref. 10) 


D @ 

D , 

*i 


9.46 lb/ft 
24.26 lb/ft 
0.61 (lb-sec )/ft 
2.09 (lb-sec)/ft 

0.72 

0.28 


(D-7) 


9: ^ 




F Y p ( D l* + K l)/ a Y x 

(D,. + K,)(W @ (e) + W c ) 

A © W C 

(D,s + K,) + A®, 

(D,s + K|) + Agj 

(D|8.K|)+A(S) 


< 8 Y t 

960082 


Figure D-2. Representation of equation (D-6). 


The equivalent weight for the controller is 

W c = 1.25 lb (D-8) 

The dynamic response of the controller was 
approximated by the spring constant, K c = 960 lb/ft. 
The measured dynamic characteristics are presented in 
appendix C. 
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APPENDIX E 

DUAL LAG EQUIVALENT TRANSFER 
FUNCTION 


frequency near 16.4 rad/sec. For one-half cycle, where 
{CO + A CO) < -60 deg/sec, the output/input would be 
as follows: 


All F-16 aircraft have a unique feedback system 
device in the forward command path referred to as a 
dual lag. This device provides for a softer roll-in and a 
sharper rollout command, thereby minimizing the ten- 
dency to overshoot. Figure E-l shows the analog sys- 
tem as excerpted from the functional block diagram 14 
describing the roll system of the F-16 XL airplane. 

A bias of ±2 0 deg/sec exists on the feedback C. 
Hence, any output within that range would be derived 
simply by computing CO = 10 RI/{s + 10). For out- 
puts greater than ±20 deg/sec, the output would be 
modified by the feedback functions to the extent of the 
output exceeding the bias. 


© ® 



Figure E-l. The nonlinear dual-lag system. 


rggll cycle = 1Q( * + 2Q) 

Ltf/J2 y (* + 3 )(j + 67) 


(E-l) 


Because the diode CD is downstream of the washout, 
for the rest of the cycle and for {CO + A CO) > -60 
deg/sec, the response would be as follows: 


'COli 

.RIP 


cycle = 


10 


{s+ 10) 


(E-2) 


Consequently, for sinusoidal inputs, a piecewise 
description of the output would be required. However, 
to avoid employing the complexity of using describing 
functions, which would result in a derivation as a func- 
tion of frequency and amplitude, a linear compromise 
was considered sufficient for determining the ampli- 
tude and phase. Figure E-2 shows this compromise 
with just the asymptotes depicting the transfer 
functions. 

The two dashed lines, as indicated, represent the two 
transfer functions 10/(5+10) and 1 0 ( v + 20)/(.y + 3) 
{s + 67). The solid line represents the linear 
compromise: 


CO = 10(5 + 40 ) 

RI ~ (5 + 6)(5 + 67 ) 


(E-3) 


The commanded roll rate and bias become important 
considerations as noted with reference to the roll ratch- 
eting experienced, described, and analyzed in this 
report (fig. 5 and 6). The oscillation is initially biased at 
an offset near —60 deg/sec. In addition, the oscillation 
varies ±10 deg/sec about the offset at the ratcheting 


At the ratcheting frequency, 16.4 rad/sec, this func- 
tion will give a linear amplitude ratio of 0.358 and a 
phase angle of -61.33°. The actual amplitude ratio is 
pointed out by the asterisk. The amplitude ratios of the 
two other functions are indicated by the small circles at 
the ratcheting frequency. 
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